Partial Differential Equations 



If the subject of ordinary differential equations is large, this is enormous. I am going to examine only one 
corner of it, and will develop only one tool to handle it: Separation of Variables. Another major tool is 
the method of characteristics and I'll not go beyond mentioning the word. When I develop a technique 
to handle the heat equation or the potential equation, don't think that it stops there. The same set of 
tools will work on the Schroedinger equation in quantum mechanics and on the wave equation in its 
many incarnations. 

10.1 The Heat Equation 

The flow of heat in one dimension is described by the heat conduction equation 



-kA 



dT 

dx 



(10.1) 



where P is the power in the form of heat energy flowing toward positive x through a wall and A is the 
area of the wall, k is the wall's thermal conductivity. Put this equation into words and it says that if a 
thin slab of material has a temperature on one side different from that on the other, then heat energy 
will flow through the slab. If the temperature difference is big or the wall is thin (dT / dx is big) then 
there's a big flow. The minus sign says that the energy flows from hot toward cold. 

When more heat comes into a region than leaves it, the temperature there will rise. This is 
described by the specific heat, c. 



dQ = mc dT, or 



dQ 
dt 



mc 



dT 
dt 



10.2) 



Again in words, the temperature rise in a chunk of material is proportional to the amount of heat added 
to it and inversely proportional to its mass. 



A 



P(x,t) 



x 



*- P(x + Ax, t) 



x + Ax 



For a slab of area A, thickness Ax, and mass density p, let the coordinates of the two sides be 
x and x + Ax. 

dQ 



m = pAAx, and 



dt 



P(x, t) - P(x + Ax, t) 



The net power into this volume is the power in from one side minus the power out from the other. Put 
these three equations together. 

dQ dT . dT . dT(x, t) . dT(x + Ax, t) 

= mc-rr = pAAxc-j- = -kA — + kA K —= ^ 

dt dt dt ox ox 

If you let Ax — > here, all you get is = 0, not very helpful. Instead divide by Ax first and then take 
the limit. 

8T_ kA / dT(x + Ax,t) dT(x,t) \ 1 
dt pcA \ dx dx J Ax 
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and in the limit this is 

8T_^(PT 
dt cp dx 2 

I was a little cavalier with the notation in that I didn't specify the argument of T on the left side. You 
could say that it was (x + Ax/2,t), but in the limit everything is evaluated at (x,t) anyway. I also 
assumed that K, the thermal conductivity, is independent of x. If not, then it stays inside the derivative, 

dT 1 d ( dT\ . ^ 

«7ir (10-4) 



In Three Dimensions 

In three dimensions, this becomes 



dt cp dx \ dx 



^ = -V 2 T (10.5) 
at cp 

Roughly speaking, the temperature in a box can change because of heat flow in any of three directions. 
More precisely, the correct three dimensional equation that replaces Eq. (10.1) is 

H = -kVT (10.6) 

where H is the heat flow vector. That is the power per area in the direction of the energy transport. 
H ■ dA = dP , the power going across the area dA. The total heat flowing into a volume is 



d Q = - jdP = - jii-dA (10.7) 



where the minus sign occurs because this is the heat flow in. For a small volume AV^, you now have 
m = pAV and 

dT Ajr dT f f* , T 
mc ^u = P AVc i^ = ~ $ H-dA 



dt h dt 

Divide by AV and take the limit as AV — > 0. The right hand side is the divergence, Eq. (9.9). 
pc^- = - lim — & H-dA = -V-H = +V ■ kVT = +kV 2 T 



dt av^o 

The last step again assumes that the thermal conductivity, K, is independent of position. 
10.2 Separation of Variables 

How do you solve these equations? I'll start with the one-dimensional case and use the method of 
separation of variables. The trick is to start by looking for a solution to the equation in the form of a 
product of a function of x and a function of t. T(x, t) = f(t)g(x). I do not assume that every solution 
to the equation will look like this — that's just not true. What will happen is that I'll be able to express 
every solution as a sum of such factored forms. That this is so is a theorem that I don't plan to prove 
here. For that you should go to a purely mathematical text on PDEs. 



10 — Partial Differential Equations 3 

If you want to find out if you have a solution, plug in: 

dT _ k d 2 T js df k <Pg 

dt cp dx 2 ' S dt® cp dx 2 

Denote the constant by n/cp = D and divide by the product fg. 

ldf _ ld 2 g 

JTt- D -gdx^ (1 °- 8) 

The left side of this equation is a function oft alone, no x. The right side is a function of x alone with 
no t, hence the name separation of variables. Because x and t can vary quite independently of each 
other, the only way that this can happen is if the two side are constant (the same constant). 

1 df „ 1 d 2 q 

- r J r = a and D--r^ = a 10.9 
j dt g dx 2 

At this point, the constant a can be anything, even complex. For a particular specified problem there 
will be boundary conditions placed on the functions, and those will constrain the a's. If a is real and 
positive then 



g{x) = AsinhJa/Dx + BcoshJa/Dx and f(t) = e at (10.10) 
For negative real a, the hyperbolic functions become circular functions. 

g(x) = A sin y/-a/D x + B cos ^J-a/D x and f(t) = e at (10.11) 

If a = then 

g(x) = Ax + B, and f(t) = constant (10.12) 

For imaginary a the f(t) is oscillating and the g(x) has both exponential and oscillatory behavior in 
space. This can really happen in very ordinary physical situations; see section 10.3. 

This analysis provides a solution to the original equation (10.3) valid for any a. A sum of such 
solutions for different a's is also a solution, for example 



T(x, t) = A ie ait sin yf-ai/D x + A 2 e a2t sin yJ-a 2 /D x 
or any other linear combination with various a's 

T(x,t) = fa{t)g a (x) 
{a's} 

It is the combined product that forms a solution to the original partial differential equation, not the 
separate factors. Determining the details of the sum is a job for Fourier series. 

Example 

A specific problem: You have a slab of material of thickness L and at a uniform temperature To. Plunge 
it into ice water at temperature T = and find the temperature inside at later times. The boundary 
condition here is that the surface temperature is zero, T(0,t) = T(L,t) = 0. This constrains the 
separated solutions, requiring that g(0) = g(L) = 0. For this to happen you can't use the hyperbolic 
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functions of x that occur when a > 0, you will need the circular functions of x, sines and cosines, 
implying that a < 0. That is also compatible with your expectation that the temperature should 
approach zero eventually, and that needs a negative exponential in time, Eq. (10.11). 



g(x) = A sin kx + B cos kx, 



with 



-a/D 



and 



/(*) 



-DkH 



g(0) = implies B = 0. g(L) = implies sin kL = 0. 

The sine vanishes for the values nn where n is any integer, positive, negative, or zero. This implies 
kL = nil, or k = nn/L. The corresponding values of a are a n = —Dn 2 ir 2 /L 2 , and the separated 
solution is 

sin (nnx/L) e -^ 2 Dt/L* (IQ.IZ) 

If n = this whole thing vanishes, so it's not much of a solution. (Not so fast there! See problem 10.2.) 
Notice that the sine is an odd function so when n < this expression just reproduces the positive n 
solution except for an overall factor of (—1), and that factor was arbitrary anyway. The negative n 
solutions are redundant, so ignore them. 

The general solution is a sum of separated solutions, see problem 10.3. 



T(*,t) = J> n sin = e-™/ L2 



(10.14) 



The problem now is to determine the coefficients a n . This is why Fourier series were invented. (Yes, 
literally, the problem of heat conduction is where Fourier series started.) At time t = you know the 
temperature distribution is T = To, a constant on < x < L. This general sum must equal To at 
time t = 0. 

oo 

1 (x, 0) = > 0"n sin —j— (0 < X < L) 
— J Lj 
1 

Multiply by sin (mnx/L^ and integrate over the domain to isolate the single term, n = m. 

-L i-L 



. mux 
ax ./ o sin — ? — = a r 



T [l — cosm7r 



L 
L 

mir 



9 mux 
ax sin — = — 



This expression for a m vanishes for even m, and when you assemble the whole series for the temperature 
you have 

///'i-/ _ ~- / ,-i i ~ 15) 



T {.,:!) 1 /;, Y ' sin ' 



m odd 



For small time, this converges, but very slowly. For large time, the convergence is very fast, often 
needing only one or two terms. As the time approaches infinity, the interior temperature approaches 
the surface temperature of zero. The graph shows the temperature profile at a sequence of times. 

T 
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The curves show the temperature dropping very quickly for points near the surface [x = or L). It 
drops more gradually near the center but eventually goes to zero. 

You can see that the boundary conditions on the temperature led to these specific boundary 
conditions on the sines and cosines. This is exactly what happened in the general development of 
Fourier series when the fundamental relationship, Eq. (5.15), required certain boundary conditions in 
order to get the orthogonality of the solutions of the harmonic oscillator differential equation. That the 
function vanishes at the boundaries was one of the possible ways to insure orthogonality. 

10.3 Oscillating Temperatures 

Take a very thick slab of material and assume that the temperature on one side of it is oscillating. Let 
the material occupy the space < x < oo and at the coordinate x = the temperature is varying in 
time as T\ cos cut. Is there any real situation in which this happens? Yes, the surface temperature of the 
Earth varies periodically from summer to winter (even in Florida). What happens to the temperature 
underground? 

The differential equation for the temperature is still Eq. (10.3), and assume that the temperature 
inside the material approaches T = far away from the surface. Separation of variables is the same 
as before, Eq. (10.9), but this time you know the time dependence at the surface. It's typical in cases 
involving oscillations that it is easier to work with complex exponentials than it is to work with sines 
and cosines. For that reason, specify that the surface temperature is T±e~ lujt instead of a cosine, 
understanding that at the end of the problem you must take the real part of the result and throw away 
the imaginary part. The imaginary part corresponds to solving the problem for a surface temperature 
of sin cut instead of cosine. It's easier to solve the two problems together then either one separately. 
(The minus sign in the exponent of e~ tujt is arbitrary; you could use a plus instead.) 

The equation (10.9) says that the time dependence that I expect is 

f dt e~ tujt K ' 

The equation for the x-dependence is then 



dx* =ag = ~ lUJ9 



This is again a simple exponential solution, say e@ x . Substitute and you have 

D/3 2 eP x = —ioje@ x , implying f3 = ±yj-iuj/D (10.16) 



Evaluate this as 

Let (3q = \JloJ2D, then the solution for the x-dependence is 

g{x) = Ae {l - i] ^ x + B e -(i-*)Aj* (10.17) 

Look at the behavior of these two terms. The first has a factor that goes as e +x and the second goes 
as e~ x . The temperature at large distances is supposed to approach zero, so that says that A = 0. 
The solutions for the temperature is now 



g e -iut e -(l-i)Pox 



(10.18) 
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The further condition is that at x = the temperature is T±e , so that tells you that B = T 



rp/ x j\ = rp ie -iut e -(l-i)P x _ rp^ e -(3 x e i(-ujt+f3 x) 
When you remember that I'm solving for the real part of this solution, the final result is 



(10.19) 




(10.20) 



This has the appearance of a temperature wave moving into the material, albeit a very strongly 
damped one. In a half wavelength of this wave, /3qX = ir, and at that point the amplitude coming from 
the exponential factor out in front is down by a factor of e~ n = 0.04. That's barely noticeable. This 
is why wine cellars are cellars. Also, you can see that at a distance where (3qX > tt/2 the temperature 
change is reversed from the value at the surface. Some distance underground, summer and winter are 
reversed. This same sort of equation comes up with the study of eddy currents in electromagnetism, 
so the same sort of results obtain. 

10.4 Spatial Temperature Distributions 

The governing equation is Eq. (10.5). For an example of a problem that falls under this heading, take a 
cube that is heated on one side and cooled on the other five sides. What is the temperature distribution 
within the cube? How does it vary in time? 

I'll take a simpler version of this problem to start with. First, I'll work in two dimensions instead 
of three; make it a very long rectangular shaped rod, extending in the ^-direction. Second, I'll look for 
the equilibrium solution, for which the time derivative is zero. These restrictions reduce the equation 
(10.5) to 

72 ^ d 2 T . d 2 T 



V 2 T 



dx 2 dy 2 







(10.21) 



I'll specify the temperature T(x,y) on the surface of the rod to be zero on three faces and To on the 
fourth. Place the coordinates so that the length of the rod is along the z-axis and the origin is in one 
corner of the rectangle. 



T(0,y) 
T(a,y) 








(0 < y < b), 
(0<y< b), 



T(x,0) 
T(x,b) 



(0 < x < a) 
T (0 < x < a) 



10.22) 



O 



x 







Look at this problem from several different angles, tear it apart, look at a lot of special cases, and 
see what can go wrong. In the process you'll see different techniques and especially a lot of applications 
of Fourier series. This single problem will illustrate many of the methods used to understand boundary 
value problems. 
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Use the same method used before for heat flow in one dimension: separation of variables. Assume 
a solution to be the product of a function of x and a function of y, then plug into the equation. 

T(x,y) = f(x)g(y), then V 2 T = < ^&-g{v) + = 

Just as in Eq. (10.8), when you divide by fg the resulting equation is separated into a term involving 
x only and one involving y only. 

l d 2 f(x) | l d 2 g(y) _ Q 
f dx 2 g dy 2 

Because x and y can be varied independently, these must be constants adding to zero. 

^d 2 f{x) , ld 2 g(y) 

f d.r 2 " and -g~W = - a (1 °- 23) 

As before, the separation constant can be any real or complex number until you start applying boundary 
conditions. You recognize that the solutions to these equations can be sines or cosines or exponentials 
or hyperbolic functions or linear functions, depending on what a is. 

The boundary conditions state that the surface temperature is held at zero on the surfaces x = 
and x = a. This suggests looking for solutions that vanish there, and that in turn says you should 
work with sines of x. In the other direction the surface temperature vanishes on only one side so you 
don't need sines in that case. The a = case gives linear functions is x and in y, and the fact that 
the temperature vanishes on x = and x = a kills these terms. (It does doesn't it?) Pick a to be a 
negative real number: call it a = —k 2 . 

- = -k 2 f =>- f (x) = Asinkx + Bcoskx 
dx 2 

The accompanying equation for g is now 

d = +k 2 g => g(y) = C sinh ky + D cosh ky 
ay 

(Or exponentials if you prefer.) The combined, separated solution to V 2 T = is 

(A sin kx + B cos kx)(C sinh ky + D cosh ky) (10.24) 

The general solution will be a sum of these, summed over various values of k. This is where you have 
to apply the boundary conditions to determine the allowed k's. 

left: T(Q,y) = = B(C sinh ky + D cosh ky), so B = 

(This holds for all y in < y < b, so the second factor can't vanish unless both C and D vanish. If 
that is the case then everything vanishes.) 

right: T(a,y) = = Asm ka(C sinh ky + D cosh ky), so sinfca = 

(The factor with y can't vanish or everything vanishes. If A = then everything vanishes. All that's 
left is sin /ca = 0.) 

bottom: T(x, 0) = = A sin kx D, so D = 
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(If A = everything is zero, so it has to be D.) 

You can now write a general solution that satisfies three of the four boundary conditions. Combine 
the coefficients A and C into one, and since it will be different for different values of k, call it 7 n . 



T ( x , V) = ^2ln sin sinh 



71=1 



nny 



(10.25) 



The mr/a appears because sinka = 0, and the limits on n omit the negative n because they are 
redundant. 

Now to find all the unknown constants 7 n , and as before that's where Fourier techniques come 
in. The fourth side, at y = b, has temperature Tq and that implies 



71=1 



. nitx . nnb 
sm sinn = 1 q 



a 



a 



On this interval < x < a these sine functions are orthogonal, so you take the scalar product of both 
side with the sine. 



f 

Jo 



dx sin 



mux 



nixx . nixb 
sm sinn 



71=1 



.%b _ r a 

a Jo 



. mnx 
dx sm 1 o 



a . rrnrb 
- 7m smh — 



To— [l-(-l)' 
mix L 



Just the odd m terms are present, m = It + 1, so the result is 

^_ 4 T y- 1 sinh ((21 + l)7ry/a) , (2f+l)7rx 
Hx, y) - n i 2^ 2e + 1 sinh , m + 1)7r6/a N sin a 



£=0 



(10.26) 



You're not done. 

Does this make sense? The dimensions are clearly correct, but after that it takes some work. 
There's really just one parameter that you have to play around with, and that's the ratio b/a. If it's 
either very big or very small you may be able to check the result. 



To 



a 0" a 

a » b b » a 

If a S> b, it looks almost like a one-dimensional problem. It is a thin slab with temperature To 
on one side and zero on the other. There's little variation along the x-direction, so the equilibrium 
equation is 

d 2 T d 2 T d 2 T 



V 2 T = 



+ 



dx 2 dy 2 dy 2 
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This simply says that the second derivative with respect to y vanishes, so the answer is the straight 
line T = A + By, and with the condition that you know the temperature at y = and at y = b you 
easily find 

T(x,y) ps T y/b 

Does the exact answer look like this? It doesn't seem to, but look closer. If b <C a then because 
< y < b you also have y <C a. The hyperbolic function factors in Eq. (10.26) will have very small 
arguments, proportional to b/a. Recall the power series expansion of the hyperbolic sine: sinhx = 
X + ■ ■ ■. These factors become approximately 

sinh ((2£ + l)iry/a) _ (2£ + l)ny/a _ y 
sinh {{2£ + l)nb/a) ~ {2£ + l)nb/a ~ b 

The temperature solution is then 

rp 1 s 4 ^ 1 y . (2£ + 1)ttx y 
T{z, y) « -T^oT+lb sm a = T °b 

1=0 

Where did that last equation come from? The coefficient of y/b is just the Fourier series of the constant 
To in terms of sines on < X < a. 

What about the opposite extreme, for which 6 > o? This is the second picture just above. 
Instead of being short and wide it is tall and narrow. For this case, look again at the arguments of the 
hyperbolic sines. Now nb/a is large and you can approximate the hyperbolic functions by going back 
to their definition. 

q x _|_ g — x \ 

sinhx = -e x , for x » 1 

2 2 ' 

The denominators in all the terms of Eq. (10.26) are large, « e^l 11 (or larger still because of the (2£ + 
1)). This will make all the terms in the series extremely small unless the numerators are correspondingly 
large. This means that the temperature stays near zero unless y is large. That makes sense. It's only 
for y near the top end that you are near to the wall with temperature To. 

You now have the case for which b^> a and y 3> a. This means that I can use the approximate 
form of the hyperbolic function for large arguments. 

S inh((2l + l)7n//q) ^ e (^y/a = ^j^,, 
sinh ((2£ + l)7r6/o) ~ e ( 2£+1 > b / a 



The temperature distribution is now approximately 

T(x, y) « ^T £ -±- e -^Mb~y)/a gin W + l)nx ^ ^ 



As compared to the previous approximation where a^> b, you can't as easily tell whether this is plausible 
or not. You can however learn from it. See also problem 10.30. 

At the very top, where y = b this reduces to the constant To that you're supposed to have at 
that position. Recall the Fourier series for a constant on < x < a. 

Move down from y = b by the distance a, so that b — y = a. That's a distance from the top 
equal to the width of the rectangle. It's still rather close to the end, but look at the series for that 
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position. 



x = x = a 

y = b 



y = b — a 



T(x, b 



71 



£=0 



2i + 1 



e -(2^ik sin 



{2i+l)nx 



For £ = 0, the exponential factor is e _ir = 0.043, and for t = 1 this factor is e" 37r = 0.00008. 
This means that measured from the To end, within the very short distance equal to the width, the 
temperature has dropped 95% of the way down to its limiting value of zero. The temperature in the 
rod is quite uniform until you are very close to the heated end. 

The Heat Flow into the Box 

All the preceding analysis and discussion was intended to make this problem and its solution sound 
oh-so-plausible. There's more, and it isn't pretty. 

The temperature on one of the four sides was given as different from the temperatures on the 
other three sides. What will the heat flow into the region be? That is, what power must you supply to 
maintain the temperature To on the single wall? 

At the beginning of this chapter, Eq. (10.1), you have the equation for the power through an 
area A, but that equation assumed that the temperature gradient dT/dx is the same all over the area 
A. If it isn't, you simply turn it into a density. 

. dT , , AP dP dT . . 

AP = ~ kAA ^ andthen AA^dA=- K dx~ (1 °- 28) 

Equivalently, just use the vector form from Eq. (10.6), H = —kVT. In Eq. (10.22) the temperature is 
To along y = b, and the power density (energy / (time ■ area)) flowing in the +y direction is —ndT/ dy, 
so the power density flowing into this area has the reversed sign, 

+ndT/dy (10.29) 

The total power flow is the integral of this over the area of the top face. 

Let L be the length of this long rectangular rod, its extent in the ^-direction. The element of 
area along the surface at y = b is then dA = Ldx, and the power flow into this face is 

a q T 

Ldx k -tt— 
o dy y=b 

The temperature function is the solution Eq. (10.26), so differentiate that equation with respect to y. 

Jo n J^o 2i+l smh((2£+l)7rb/a) a y 
/ dx > i 



4LkT f a , ^ (2£ + l)nx 
i ./.. \ gm 



£=0 
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and this sum does not converge. I'm going to push ahead anyway, temporarily pretending that I didn't 
notice this minor difficulty with the series. Just go ahead and integrate the series term by term and 
hope for the best. 



This infinite series for the total power entering the top face is infinite. The series doesn't converge (use 
the integral test). 

This innocuous-seeming problem is suddenly pathological because it would take an infinite power 
source to maintain this temperature difference. Why should that be? Look at the corners. You're 
trying to maintain a non-zero temperature difference (To — 0) between two walls that are touching. 
This can't happen, and the equations are telling you so! It means that the boundary conditions specified 
in Eq. (10.22) are impossible to maintain. The temperature on the boundary at y = b can't be constant 
all the way to the edge. It must drop off to zero as it approaches x = and x = a. This makes the 
problem more difficult, but then reality is typically more complicated than our simple, idealized models. 

Does this make the solution Eq. (10.26) valueless? No, it simply means that you can't push it 
too hard. This solution will be good until you get near the corners, where you can't possibly maintain 
the constant-temperature boundary condition. In other regions it will be a good approximation to the 
physical problem. 

10.5 Specified Heat Flow 

In the previous examples, I specified the temperature on the boundaries and from that I determined the 
temperature inside. In the particular example, the solution was not physically plausible all the way to 
the edge, though the mathematics were (I hope) enlightening. Instead, I'll reverse the process and try 
to specify the size of the heat flow, computing the resulting temperature from that. This time perhaps 
the results will be a better reflection of reality. 

Equation (10.29) tells you the power density at the surface, and I'll examine the case for which 
this is a constant. Call it Fq. (There's not a conventional symbol, so this will do.) The plus sign occurs 
because the flow is into the box. 



The other three walls have the same zero temperature conditions as Eq. (10.22). Which forms of the 
separated solutions must I use now? The same ones as before or different ones? 
Look again at the a = solutions to Eqs. (10.23). That solution is 



In order to handle the fact that the temperature is zero at y = and that the derivative with respect 
to y is given at y = b, 





+K 7 -(x,b) = F 



(A + Bx)(C + Dy) 



{A + Bx){C) 







and 



(A + Bx)(D) 



F /k, 



implying 



C = = B 



then 



AD 




(10.30) 



This matches the boundary conditions at both y = and y = b. All that's left is to make everything 
work at the other two faces. 
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• y l 1111 
















Foy/it 





If I can find a solution that equals —F^y/n on the left and right faces then it will cancel the 
+Fqu/k that Eq. (10.30) provides. But I can't disturb the top and bottom boundary conditions. The 
way to do that is to find functions that equal zero at y = and whose derivative equals zero at y = b. 
This is a familiar sort of condition that showed up several times in chapter five on Fourier series. It 
is equivalent to saying that the top surface is insulated so that heat can't flow through it. You then 
use superposition to combine the solution with uniform heat flow and the solution with an insulated 
boundary. 

Instead of Eq. (10.24), use the opposite sign for a, so the solutions are of the form 

(^4 sin ky + B cos ky) (C sinh kx + D cosh kx) 

I require that this equals zero at y = 0, so that says 

(0 + B) (C sinh kx + D cosh kx) = 

so B = 0. Now require that the derivative equals zero at y = b, so 

Ak cos kb = 0, or kb = (n + \)n for n = 0, 1, 2 . . . 

The value of the temperature is the same on the left that it is on the right, so 

C sinh kO + D cosh kO = C sinh ka + D cosh ka =^ C = D(l — cosh ka)/ sinh ka (10.31) 

This is starting to get messy, so it's time to look around and see if I've missed anything that 
could simplify the calculation. There's no guarantee that there is any simpler way, but it is always worth 
looking. The fact that the system is the same on the left as on the right means that the temperature 
will be symmetric about the central axis of the box, about x = a/2. That it is even about this point 
implies that the hyperbolic functions of x should be even about x = a/2. You can do this simply by 
using a cosh about that point. 

A sin ky(D cosh k(x - §)) 
Put these together and you have a sum 

(n+ \)ny 



n=0 



sin 



cosh 



(n + \)tx{x - 1) 



(10.32) 



Each of these terms satisfies Laplace's equation, satisfies the boundary conditions at y = and y = b, 
and is even about the centerline x = a/2. It is now a problem in Fourier series to match the conditions 
at x = 0. They're then automatically satisfied at x = a. 



n=0 



a n sin 



(n+\)iry 



cosh 



(n + \)-Ka 



2b 



f y 

K 



(10.33) 
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The sines are orthogonal by the theorem Eq. (5.15), so you can pick out the component a n by the 
orthogonality of these basis functions. 



u, 



sm 



(n+ \)ixy 



then 



or, a m (u m ,Um) cosh 



[u m , left side) = (u m , right side) 



(m + \)na 
2b 



F 



(u m ,y) 



Write this out; do the integrals, add the linear term, and you have 

Mr „ 



T{x,y) = F --^^ 



sm 



(2n + l) 2 



n=0 

(n + |)7ry N 



cosh 



(n + i)7r(x 



(10.34) 



scch 



(n + |)7ra 
26 



Now analyze this to see if it makes sense. I'll look at the same cases as the last time: b <C a 
and a <C 6. The simpler case, where the box is short and wide, has b <^ a. This makes the arguments 
of the cosh and sech large, with an a/b in them. For large argument you can approximate the cosh by 

cosh x « e x /2, x » 1 

Now examine a typical term in the sum (10.34), and I have to be a little more specific and choose x on 
the left or right of a/2. The reason for that is the preceding equation requires x large and positive. I'll 
take x on the right, as it makes no difference. The hyperbolic functions in (10.34) are approximately 



exp ((n + 1)tt(x - §)/&) 
exp ((n + \)na/2b) 



,((2n+l)7r(a:-a)/26) 



As long as X is not near the end, that is, not near x = a, the quantity in the exponential is large and 
negative for all n. The exponential in turn makes this extremely small so that the entire sum becomes 
negligible. The temperature distribution is then the single term 



T(x,y)^F Q 



It's essentially a one dimensional problem, with the heat flow only along the — y direction. 

In the reverse case for which the box is tall and thin, a <C b, the arguments of the hyperbolic 
functions are small. This invites a power series expansion, but that approach doesn't work. The analysis 
of this case is quite tricky, and I finally concluded that it's not worth the trouble to write it up. It leads 
to a rather complicated integral. 

10.6 Electrostatics 

The equation for the electrostatic potential in a vacuum is exactly the same as Eq. (10.21) for the 
temperature in static equilibrium, V 2 V = 0, with the electric field E = — VV. The same equation 
applies to the gravitational potential, Eq. (9.42). 

Perhaps you've looked into a microwave oven. You can see inside it, but the microwaves aren't 
supposed to get out. How can this be? Light is just another form of electromagnetic radiation, so why 



10 — Partial Differential Equations 



14 



does one EM wave get through while the other one doesn't? I won't solve the whole electromagnetic 
radiation problem here, but I'll look at the static analog to get some general idea of what's happening. 



///// 



-L L 2L x 

Arrange a set of conducting strips in the x-y plane and with insulation between them so that 
they don't quite touch each other. Now apply voltage Vo on every other one so that the potentials are 
alternately zero and Vq. This sets the potential in the z = plane to be independent of y and 

z = 0: V (x,y) = {V° |°<x<2L) V(x + 2L, y) = V(x, y), all x, y (10.35) 

What is then the potential above the plane, z > 0? Above the plane V satisfies Laplace's equation, 

„ 2l , d 2 V d 2 V d 2 V 

W = -3-5- + -3-5- + -3-5- = 10.36 
ox 2 ay 1 oz 2 

The potential is independent of y in the plane, so it will be independent of y everywhere. Separate 
variables in the remaining coordinates. 

This is separated as a function of x plus a function of y, so the terms are constants. 

I've chosen the separation constant in this form because the boundary condition is periodic in x, and 
that implies that I'll want oscillating functions there, not exponentials. 

f( x ) = e iax and f( x + 2L) = f{x) 
e 2Lia = 1, or 2La = 2nn, n = 0, ±1, ±2,... 

The separated solutions are then 

f(x)g(z) = e ™ix/L^ Ae n-Kz/L + Be ^z/L^ (10.38) 

The solution for z > is therefore the sum 

00 

V(x,z)= e nnix / L (A n e nnz / L + B n e' nnz / L ) (10.39) 

n=— 00 

The coefficients A n and B n are to be determined by Fourier techniques. First however, look at the 
z-behavior. As you move away from the plane toward positive z, the potential should not increase 
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without bound. Terms such as e 712 /^ however increase with z. This means that the coefficients of the 
terms that increase exponentially in z cannot be there. 

A n = for n > 0, and B n = for n < 

oo —1 
V(x, Z) = A + B + e nnix / L B n e~ n7rz / L + Y ^ix/L^nnz/L ( 1Q 40 ) 

n=l n=— oo 

The combined constant Aq + Bq is really one constant; you can call it Co if you like. Now use the 
usual Fourier techniques given that you know the potential at z = 0. 

oo —1 

V(x, 0) = Co + B n e n7rix/L + Y ^ne nnix/L 

n=l n=— oo 

The scalar product of e miTix /L w j tn this equation is 

f 2LC (m = 0) 
(e m7vix / L , V(x, 0)> = I 2LB m (m > 0) 
{ 2LA m (m < 0) 

Now evaluate the integral on the left side. First, m^0: 



(10.41) 



r-L 

(e mnix / L ,V(x,0)) =j dxe 



c/L /O (-L < x < 0) 
V (0 < x < L) 



V I dxe~ mnix / L = V — —. e 
Jo -mm 



-rmrix/L 



V - 



-VTLTXl 



j[(-l) m -l] 



Then evaluate it separately for m = 0, and you have (l, V(x, 0)) = VqL. 

Now assemble the result. Before plunging in, look at what will happen. 
The m = term sits by itself. 

For the other terms, only odd m have non-zero values. 

V(x, z) = -V + V J2 [(-!) m - 1] e mnix / L e- mnz / L 

Lit I i l\ i 



m=l 

-1 



(10.42) 



+ Vb I rr 1 \m_-n e mmx/L e +m-Kz/L 

^ -2rrmi 1 J 

m=-oo 

To put this into a real form that is easier to interpret, change variables, letting m = —n in the second 
sum and m = n in the first, finally changing the sum so that it is over just the odd terms. 

1 00 1 

V(x, z) = l -V + V £ [{-XT - 1] e n ™' L e-<"«/ L 

n=l 

oo 

+ y V- r/^n _ -n e -n7rix/L e -nnz/L 

^ +2nm 1 J 

1 00 1 

= -Vo + V Y \(-l) n - ll — sm(n7rx/L)e- n7VZ / L 

n=l 

1 9 °° 1 



(10.43) 
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Having done all the work to get to the answer, what can I learn from it? 
What does it look like? 

Are there any properties of the solution that are unexpected? 
Should I have anticipated the form of the result? 
Is there an easier way to get to the result? 

To see what it looks like, examine some values of z, the distance above the surface. If z = L, 
the coefficient for successive terms is 

I = : -e- 71 = 0.028 £ = 1 : — = 1.7 x 10~ 5 (10.44) 

7T 67T 

The constant term is the average potential, and the £ = term adds only a modest ripple, about 5% of 
the constant average value. If you move up to z = 2L the first factor is 0.0012 and that's a little more 
than 0.2% ripple. The sharp jumps from +Vo to zero and back disappear rapidly. That the oscillations 
vanish so quickly with distance is perhaps not what you would guess until you have analyzed such a 
problem. 



■V(x,L/2) 

*-V(x,0) 

The graph shows the potential function at the surface, z = 0, as it oscillates between Vq and 
zero. It then shows successive graphs of Eq. (10.43) at z = L/2, then at z = L, then at z = 1.5L. 
The ripple is barely visible at the that last distance. The radiation through the screen of a microwave 
oven is filtered in much the same way because the wavelength of the radiation is large compared to the 
size of the holes in the screen. 

When you write the form of the series for the potential, Eq. (10.40), you can see this coming if 
you look for it. The oscillating terms in x are accompanied by exponential terms in z, and the rapid 
damping with distance is already apparent: e - n7VZ /L^ y ou don't have to solve for a single coefficient 
to see that the oscillations vanish very rapidly with distance. 

The original potential on the surface was neither even nor odd, but except for the constant 
average value, it is an odd function. 



z = : V(x, y) 



l V ° + {-vJ/l £<»< % V {x + 2L,y ) = V {x ,y) (10.45) 



Solve the potential problem for the constant Vq/2 and you have a constant. Solve it for the remaining 
odd function on the boundary and you should expect an odd function for V(x,z). If you make these 
observations before solving the problem you can save yourself some algebra, as it will lead you to the 
form of the solution faster. 

The potential is periodic on the x-y plane, so periodic boundary conditions are the appropriate 
ones. You can express these in more than one way, taking as a basis for the expansion either complex 
exponentials or sines and cosines. 

nnix/L „ _ n + -i 

e n-u,±i,... (1Q46) 

or the combination cos(nnx/L), n = 0, 1, . . . sm(niix/L), n = 1, 2, . . . 

For a random problem with no special symmetry the exponential choice typically leads to easier integrals. 
In this case the boundary condition has some symmetry that you can take advantage of: it's almost 
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odd. The constant term in Eq. (10.30) is the n = element of the cosine set, and that's necessarily 
orthogonal to all the sines. For the rest, you do the expansion 

-Vq/2 (0<x<L) A s Wmrr/n 
-Vo/2 (L<x<2L) -Z^^sm(n7rx/L) 

The odd term in the boundary condition (10.45) is necessarily a sum of sines, with no cosines. The 
cosines are orthogonal to an odd function. See problem 10.11. 

More Electrostatic Examples 

Specify the electric potential in the x-y plane to be an array, periodic in both the x and the y-directions. 
V(x,y,z = 0) is Vo on the rectangle (0 < X < a, < y < b) as well as in the darkened boxes in the 
picture; it is zero in the white boxes. What is the potential for z > 0? 




The equation is still Eq. (10.36), but now you have to do the separation of variables along all 
three coordinates, V(x,y,z) = f(x)g(y)h(z). Substitute into the Laplace equation and divide by 

fgh. 

1 d 2 f 1 d 2 g 1 d 2 h 



+ 







/ dx 2 g dy 2 h dz 2 

These terms are functions of the single variables x, y, and z respectively, so the only way this can work 
is if they are separately constant. 



fdx 2 



k 2 



ld?g 
gdy 2 



-k 2 



2- 



1 d 2 h _ 2 2 _ 2 
hJz^~ kx+k2 ~^ 



I made the choice of the signs for these constants because the boundary function is periodic in x and 
in y, so I expect sines and cosines along those directions. The separated solution is 



(A sin k\X + B cos k\x){C sin k 2 y + D cos k 2 y)(Ee k3Z + Fe 



-k 3 z\ 



(10.47) 



What about the case for separation constants of zero? Yes, that's needed too; the average value of 
the potential on the surface is Vo/2, so just as with the example leading to Eq. (10.43) this will have 
a constant term of that value. The periodicity in x is 2a and in y it is 2b, so this determines 



ki=nir/a, k2=m7r/b 



then 



2tt2 



1, 2,... 



where n and m are independent integers. Use the experience that led to Eq. (10.45) to write V on the 
surface as a sum of the constant Vo/2 and a function that is odd in both x and in y. As there, the 
odd function in x will be represented by a sum of sines in x, and the same statement will hold for the 
y coordinate. This leads to the form of the sum 



V(x,y,z) = ^Vo + Y^ 

n=l m=l 



,,,„ ( — Jsm( - . - ]<■-' 
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where k nm is the of the preceding equation. What happened to the other term in z, the one with 
the positive exponent? Did I say that I'm looking for solutions in the domain z > 0? 

At z = this must match the boundary conditions stated, and as before, the orthogonality of the 
sines on the two domains allows you to determine the coefficients. You simply have to do two integrals 
instead of one. See problem 10.19. 

odd n odd m 



10.7 Cylindrical Coordinates 

Rectangular coordinates are not always the right choice. Cylindrical, spherical, and other choices are 
often needed. For cylindrical coordinates, the gradient and divergence are, from Eqs. (9.24) and (9.15) 

„ T , „dV 2 ldV „dV , „ _ ld{rv r ) ldv^ dv z 

or r o<p oz r or r o<p dz 

Combine these to get the Laplacian in cylindrical coordinates. 

7 2 Tr 1 d ( dV\ . 1 d 2 V . d 2 V 



S/.S/V = V 2 V = -j-{r^-\+^— 7I + — T (10.49) 
r Or \ Or J r z o<p A oz 1 

For electrostatics, the equation remains V 2 V = 0, and you can approach it the same way as before, 
using separation of variables. I'll start with the special case for which everything is independent of z. 
Assume a solution of the form V = f(r)g((f>), then 

r dr \ dr ) r 2 d<p 2 
Multiply this by r 2 and divide by f{r)g(<j)) to get 

r d ( df\ 1 d 2 g 



f dr \ dr J g d(p 2 

This is separated. The first terms depends on r alone, and the second term on alone. For this to 
hold the terms must be constants. 



rdfdf\_ and ld?g 
f dr \ dr J g d<p 



rj-)=a and = -a (10.50) 



The second equation, in 0, is familiar. If a is positive this is a harmonic oscillator, and that is the most 
common way this solution is applied. I'll then look at the case for a > 0, for which the substitution 
a = n 2 makes sense. 

a = 0: = g(<j>) = A + B<f) 

4r(4)=°^ m = C + D*r 
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d?q 

a = n 2 > = ~~ n2 9 ==> * g(4>) = A cos n<J) + B sin n<J) 

r 2 ^+r^-n 2 f = => f(r) = Cr n + D r - n 

There's not yet a restriction that n is an integer, though that's often the case. Verifying the last 
solution for / is easy. 

A general solution that is the sum of all these terms is 

V(r,4>) = (C Q + D o \nr)(A o + B o( p) + ^(C n r n + D n r- n )(A nC osn4> + B n smn0) (10.51) 

n 

Some of these terms should be familiar: 
CqAq is just a constant potential. 

-Do^4ol nr is the potential of a uniform line charge; dlnr/dr = l/r, and that is the way that the 
electric field drops off with distance from the axis. 

CiAircos^ is the potential of a uniform field (as is the rsin0 term). Write this in the form 
C\A\T cos0 = C\A\X, and the gradient of this is C\A\X. The sine gives y. 
See problem 10.24. 

Example 

A conducting wire, radius R, is placed in a uniform electric field Eq, and perpendicular to it. Put the 
wire along the z-axis and call the positive x-axis the direction that the field points. That's = 0. 
In the absence of the wire, the potential for the uniform field is V = —EqX = — E^r cos (ft, because 
— VV = EqX. The total solution will be in the form of Eq. (10.51). 

Now turn the general form of the solution into a particular one for this problem. The entire range 
of from to 2tc appears here; you can go all the way around the origin and come back to where you 
started. The potential is a function, meaning that it's single valued, and that eliminates Bq(J). It also 
implies that all the n are integers. The applied field has a potential that is even in 0. That means that 
you should expect the solution to be even in too. Does it really? You also have to note that the 
physical system and its attendant boundary conditions are even in — it's a cylinder. Then too, the 
first few times that you do this sort of problem you should see what happens to the odd terms; what 
makes them vanish? I won't eliminate the sin0 terms from the start, but I'll calculate them and show 
that they do vanish. 

oo 

V(r, 0) = (Co + A) In r ) B + J] (C n r n + D n r- n )(A n cos n<P + B n sin mj>) 

71=1 

Carrying along all these products of (still unknown) factors such as D n A n is awkward. It makes it look 
neater and it is easier to follow if I combine and rename some of them. 

oo oo 

V(r, 0) = C + D \nr + ^(C n r n + D n r~ n ) cosn0 + ^(C^r 71 + D' n r~ n ) sinn0 (10.52) 

n=l n=l 

As r — > oo, the potential looks like —E$r cos 0. That implies that C n = for n > 1, and that C' n = 
for all n, and that C\ = —Eq. 

Now use the boundary conditions on the cylinder. It is a conductor, so in this static case the 
potential is a constant all through it, in particular on the surface. I may as well take that constant to 
be zero, though it doesn't really matter. 

oo oo 

V(R, 0) = = C + D In R + ^ {C n R n + D n R~ n ) cos n<f) + ^ (C' n R n + D' n R- n ) sin n<f) 

n=l n=l 
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Multiply by sinm0 and integrate over 0. The trigonometric functions are orthogonal, so all that 
survives is 

= {C' m R m + D' m R- m )7r all m > 1 

That gets rid of all the rest of the sine terms as promised: D' m = for all m because C' m = for all 
m. Now repeat for cosm0. 

= C + D In R (m = 0) and = (C m R m + D m R" m )n (m > 0) 

All of the C m = for m > 1, so this says that the same applies to D m . The m = 1 equation 
determines L>i in terms of Ci and then Eq. 

D x = -C X R 2 = +E R 2 

Only the Co and Dq terms are left, and that requires another boundary condition. When specifying 
the problem initially, I didn't say whether or not there is any charge on the wire. In such a case you 
could naturally assume that it is zero, but you have to say so explicitly because that affects the final 
result. Make it zero. That kills the lnr term. The reason for that goes back to the interpretation of 
this term. Its negative gradient is the electric field, and that would be — l/r, the field of a uniform line 
charge. If I assume there isn't one, then Dq = and so the same for Co. Put this all together and 

R 2 

V{r,<j>) = -E r cos (p + E — cos0 (10.53) 
The electric field that this presents is, from Eq. (9.24) 

E = —VV = E (f cos <p — 4> sin <p) — EqR 2 ^— f cos <p — 0^ sin < 

R 2 

= EqX + Eq—^{t cos0 + 0sin0) 
As a check to see what this looks like, what is the electric field at the surface of the cylinder? 

E(R, 0) = Eo(f cos <fi — 0sin0) — EqR 2 ^— f-^ cos0 — sin0^ = 2Eofcos<p 



It's perpendicular to the surface, as it should be. At the left and right, = 0, 7T, it is twice as large as 
the uniform field alone would be at those points. 
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Problems 

10.1 The specific heat of a particular type of stainless steel (CF8M) is 500J/kg-K. Its thermal 
conductivity is 16.2 W/m-K and its density is 7750 kg/m 3 . A slab of this steel 1.00 cm thick is at 
a temperature 100°C and it is placed into ice water. Assume the simplest boundary condition that its 
surface temperature stays at zero, and find the internal temperature at later times. When is the 2 nd 
term in the series, Eq. (10.15), only 5% of the 1 st ? Sketch the temperature distribution then, indicating 
the scale correctly. 

10.2 In Eq. (10.13) I eliminated the n = solution by a fallacious argument. What is a in this case? 
This gives one more term in the sum, Eq. (10.14). Show that with the boundary conditions stated, this 
extra term is zero anyway (this time). 

10.3 In Eq. (10.14) you have the sum of many terms. Does it still satisfy the original differential 
equation, Eq. (10.3)? 

10.4 In the example Eq. (10.15) the final temperature was zero. What if the final temperature is T\l 
Or what if I use the Kelvin scale, so that the final temperature is 273°? Add the appropriate extra term, 
making sure that you still have a solution to the original differential equation and that the boundary 
conditions are satisfied. 

10.5 In the example Eq. (10.15) the final temperature was zero on both sides. What if it's zero on 
just the side at x = while the side at x = L stays at To? What is the solution now? 

Ans: Tqx/L + (2T /tt) £f (l/ra) sm(mrx/L)e~ n2lv2Dt / L2 

10.6 You have a slab of material of thickness L and at a uniform temperature To. The side at x = L is 
insulated so that heat can't flow in or out of that surface. By Eq. (10.1), this tells you that dT / dx = 
at that surface. Plunge the other side into ice water at temperature T = and find the temperature 
inside at later time. The boundary condition on the x = surface is the same as in the example in the 
text, T(0, t) = 0. (a) Separate variables and find the appropriate separated solutions for these boundary 
conditions. Are the separated solutions orthogonal? Use the techniques of Eq. (5.15). (b) When the 
lowest order term has dropped to where its contribution to the temperature at x = L is To/2, how big 
is the next term in the series? Sketch the temperature distribution in the slab at that time. 

Ans: (4T / 7 r)X; "(2^p T )sm [(n + \)t:x/ L]e^ n+l l 2 ^ 2Dt l L \ -9.43 x 10~ 5 To 

10.7 In the analysis leading to Eq. (10.26) the temperature at y = b was set to To. If instead, you 
have the temperature at x = a set to To with all the other sides at zero, write down the answer for 
the temperature within the rod. Now use the fact that Eq. (10.21) is linear to write down the solution 
if both the sides at y = b and x = a are set to To. 

10.8 In leading up to Eq. (10.25) I didn't examine the third possibility for the separation constant, 
that it's zero. Do so. 

10.9 Look at the boundary condition of Eq. (10.22) again. Another way to solve this problem is to use 
the solution for which the separation constant is zero, and to use it to satisfy the conditions at y = 
and y = b. You will then have one term in the separated solution that is T^y/b, and that means that 
in Eq. (10.23) you will have to choose the separation variable to be positive instead of negative. Why? 
Because now all the rest of the terms in the sum over separated solutions must vanish at y = and 
y = b. You've already satisfied the boundary conditions on those surfaces by using the T$y/b term. 
Now you have to satisfy the boundary conditions on x = and x = a because the total temperature 
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there must be zero. That in turn means that the sum over all the rest of the separated terms must add 
to —Toy/b at x = and x = a. When you analyze this solution in the same spirit as the analysis of 
Eq. (10.26), compare the convergence properties of that solution to your new one. In particular, look 
at a -C b and a S> b to see which version converges better in each case. 
Ans: T y/b+(2T /n)Y,T [{- l T/ n ] sm(nny/b) cosh (nir(x - a/2)/b)/ cosh(mra/2&) 

10.10 Finish the re-analysis of the electrostatic boundary value problem Eq. (10.45) starting from 
Eq. (10.46). This will get the potential for z ^ with perhaps less work than before. 

10.11 Examine the solution Eq. (10.42) at z = in the light of problem 5.11. 

10.12 A thick slab of material is alternately heated and cooled at its surface so the its surface tem- 
perature oscillates as 



T(0,t) 



Tx (0 < t < t ) 
-Ti {t <t<2t ) 



T(0,t + 2t )=T(0,t) 



That is, the period of the oscillation is 2to- Find the temperature inside the material, for x > 0. 

does this behavior differ from the solution in Eq. (10.20)? 

Ans: ^Er=o( 1 /(2fc + l))e- /3 ^sin((2A; + lM-/3 fc x)); u = n/t 



How 



10.13 Fill in the missing steps in finding the solution, Eq. (10.34). 

10.14 A variation on the problem of the alternating potential strips in section 10.6. Place a grounded 
conducting sheet parallel to the x-y plane at a height z = d above it. The potential there is then 
V(x, y,z = d) = 0. Solve for the potential in the gap between z = and z = d. A suggestion: you 
may find it easier to turn the coordinate system over so that the grounded sheet is at z = and the 
alternating strips are at z = d. This switch of coordinates is in no way essential, but it is a bit easier. 
Also, I want to point out that you will need to consider the case for which the separation constant in 
Eq. (10.37) is zero. 

10.15 Starting from Eq. (10.52) and repeat the example there, but assume that the conducting wire 
is in an external electric field E^y instead of EqX. Repeat the calculation for the potential and for the 
electric field, filling in the details of the missing steps. 

10.16 A very long conducting cylindrical shell of radius R is split in two along lines parallel to its axis. 
The two halves are wired to a circuit that places one half at potential Vq and the other half at potential 
— Vo- (a) What is the potential everywhere inside the cylinder? Use the results of section 10.7 and 
assume a solution of the form 



V(r, 4>) = ^2 r?1 { a n cos n + b n sin n(, 




Match the boundary condition that 



V(R,< 



V (0 < (p < it) 
-V (it <4>< 2n) 
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I picked the axis for = pointing toward the split between the cylinders. No particular reason, but 
you have to make a choice. I make the approximation that the cylinder is infinitely long so that z 
dependence doesn't enter. Also, the two halves of the cylinder almost touch so I'm neglecting the 
distance between them. 

(b) What is the electric field, — W on the central axis? Is this answer more or less what you would 
estimate before solving the problem? Ans: (b) E = ^Vq/tiR. 

10.17 Solve the preceding problem outside the cylinder. The integer n can be either positive or 
negative, and this time you'll need the negative values. (And why must n be an integer?) 

Ans: (4Vo/7r)E.odd(l/^(^/O n smn0 

10.18 In the split cylinder of problem 10.16, insert a coaxial wire of radius R\ < R. It is at zero 
potential. Now what is the potential in the domain R\ <r < Rl You will need both the positive and 
negative n values, ^(^n^ 71 + B n r~ n ) smncj) 

Ans: (4Vo/tt) Em odd sin rri(f)[-R^ m r m + R™ r - m }/m{R~ m Rf - R^ m R m ] 

10.19 Fill in the missing steps in deriving Eq. (10.48). 

10.20 Analyze how rapidly the solution Eq. (10.48) approaches a constant as z increases from zero. 
Compare Eq. (10.44). 

10.21 A broad class of second order linear homogeneous differential equations can, with some manip- 
ulation, be put into the form (Sturm-Liouville) 

(p(x)u' )' + q(x)u = \w(x)u 

Assume that the functions p, q, and w are real, and use manipulations much like those that led to 
the identity Eq. (5.15). Derive the analogous identity for this new differential equation. When you use 
separation of variables on equations involving the Laplacian you will commonly come to an ordinary 
differential equation of exactly this form. The precise details will depend on the coordinate system you 
are using as well as other aspects of the PDE. 

10.22 Carry on from Eq. (10.31) and deduce the separated solution that satisfies these boundary 
condition. Show that it is equivalent to Eq. (10.32). 

10.23 The Laplacian in cylindrical coordinates is Eq. (10.49). Separate variables for the equation 
V 2 V = and you will see that the equations in z and are familiar. The equation in the r variable 
is less so, but you've seen it (almost) in Eqs. (4.18) and (4.22). Make a change of variables in the 
r-differential equation, r = kr' , and turn it into exactly the form described there. 

10.24 In the preceding problem suppose that there's no ^-dependence. Look at the case where the 
separation constant is zero for both the r and functions, finally assembling the product of the two 
for another solution of the whole equation. 

These results provide four different solutions, a constant, a function of r alone, a function of alone, 
and a function of both. In each of these cases, assume that these functions are potentials V and that 
E = —W is the electric field from each potential. Sketch equipotentials for each case, then sketch 
the corresponding vector fields that they produce (a lot of arrows). 
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10.25 Do problem 8.23 and now solve it, finding all solutions to the wave equation. 
Ans: f(x — vt) + g(x + vt) 

10.26 Use the results of problem 10.24 to find the potential in the corner between two very 
large metal plates set at right angles. One at potential zero, the other at potential Vq. 
Compute the electric field, — VV and draw the results. Ans: —2Vo(J)/irr 

10.27 A thin metal sheet has a straight edge for one of its boundaries. Another thin metal sheet is 
cut the same way. The two straight boundaries are placed in the same plane and almost, but not quite 
touching. Now apply a potential difference between them, putting one at a voltage Vq and the other at 
—Vq. In the region of space near to the almost touching boundary, what is the electric potential? From 
that, compute and draw the electric field. 

10.28 A slab of heat conducting material lies between coordinates x = —L and x = +L, which are 
at temperatures T\ and T2 respectively. In the steady state (dT/dt = 0), what is the temperature 
distribution inside? Now express the result in cylindrical coordinates around the z-axis and show how 
it matches the sum of cylindrical coordinate solutions of V 2 T = from problem 10.15. What if the 
surfaces of the slab had been specified at y = —L and y = +L instead? 

10.29 The result of problem 10.16 has a series of terms that look like (x n /n) smncf) (odd n). You 
can use complex exponentials, do a little rearranging and factoring, and sum this series. Along the way 
you will have to figure out what the sum z + z 3 /3 + z 5 /5 + • • • is. Refer to section 2.7. Finally of 
course, the answer is real, and if you look hard you may find a simple interpretation for the result. Be 
sure you've done problem 10.24 before trying this last step. Ans: 2Vq{6\ + 62) /ft. You still have to 
decipher what 6\ and 62 are. 

10.30 Sum the series Eq. (10.27) to get a closed-form analytic expression for the temperature distri- 
bution. You will find the techniques of section 5.7 useful, but there are still a lot of steps. Recall also 

ln(r e t0 ) = lnr + id. Ans: (2Tq/tt) tan -1 [sm(nx/a)/ sinh(7r(6 — y)/a)] 

10.31 A generalization of the problem specified in Eq. (10.22). Now the four sides have temperatures 
given respectively to be the constants T\, T2, T3, T4. Note: with a little bit of foresight, you won't 
have to work very hard at all to solve this. 

10.32 Use the electrostatic equations from problem 9.21 and assume that the electric charge density is 
given by p = poa/r, where this is in cylindrical coordinates, (a) What cylindrically symmetric electric 
field comes from this charge distribution? (b) From E = — W what potential function V do you get? 

10.33 Repeat the preceding problem, but now interpret r as referring to spherical coordinates. What 

is V 2 V7 

10.34 The Laplacian in spherical coordinates is Eq. (9.43). The electrostatic potential equation is 
V 2 V = just as before, but now take the special case of azimuthal symmetry so that the potential 
function is independent of (f). Apply the method of separation of variables to find solutions of the form 
f(r)g(9). You will get two ordinary differential equations for / and g. The second of these equations 
is much simpler if you make the change of independent variable x = cos 9. Use the chain rule a couple 
of times to do so, showing that the two differential equations are 
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10.35 For the preceding equations show that there are solutions of the form f(r) = Ar n , and recall 
the analysis in section 4.11 for the g solutions. What values of the separation constant C will allow 
solutions that are finite as x — >■ ±1 (9 — > 0, 7r)? What are the corresponding functions of r? Don't 
forget that there are two solutions to the second order differential equation for / — two roots to a 
quadratic equation. 

10.36 Write out the separated solutions to the preceding problem (the ones that are are finite as 9 
approaches or 7r) for the two smallest allowed values of the separation constant C: and 2. In each 
of the four cases, interpret and sketch the potential and its corresponding electric field, — VV. How do 
you sketch a potential? Draw equipotentials. 

10.37 From the preceding problem you can have a potential, a solution of Laplace's equation, in the 
form (Ar + B/r 2 ) cos9. Show that by an appropriate choice of A and B, this has an electric field that 
for large distances from the origin looks like Eq2, and that on the sphere r = R the total potential is 
zero — a grounded, conducting sphere. What does the total electric field look like for r > R; sketch 
some field lines. Start by asking what the electric field is as r — > R. 

10.38 Similar to problem 10.16, but the potential on the cylinder is 

V(t> ja _ fVo (0 < < 7r/2 and n < </>< 3n/2) 
[ ' ^ ~ \ -V Q (tt/2 < < 7T and 3tt/2 < < 2tt) 

Draw the electric field in the region near r = 0. 

10.39 What is the area charge density on the surface of the cylinder where the potential is given by 
Eq. (10.53)? 



